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It  is  known  that  a  spatially  localized  rf  electric  field  E  of  frequency  w 
produces  an  average  static  force,  i.e.,  a  ponderomotive  force.*  In  the  pres¬ 
ence  of  a  static  magnetic  field  B  a  large  enhancement  of  the  ponderomotive 
2 

force  is  predicted  to  appear  as  w  approaches  the  particle  gyrofrequency  Q. 

2 

The  standard  expression  given  for  the  ponderomotive  (pseudo)  potential  is 

2  |  E  | 2 

*p  -  —*2 . ~T~  (1) 

where  q  and  m  are  the  charge  and  mass  of  the  particle.  The  non-physical 
singularity  in  Eq.  (1)  arises  because  in  obtaining  this  expression  the  grad¬ 
ient  scale  length  SL  of  the  electric  field  is  assumed  to  be  large  compared 
to  v/(w  -  n) ,  where  v  is  the  unperturbed  particle  velocity.  However,  this 

3 

adiabatic  approximation  is  not  valid  near  gyroresonance,  hence  the  singu¬ 
larity  in  Eq.  (1)  is  suspect.  Since  the  ponderomotive  force  plays  a  funda- 
4 

mental  role  in  many  nonlinear  processes  in  plasma  physics,  it  is  important 

to  understand  the  singular  behavior  near  the  gyrofrequency.  In  particular, 

efforts  are  in  progress  to  utilize  the  singular  nature  of  ^  to  selectively 

confine  ions  of  different  q/m  for  isotope  separation^  6  and  to  plug  open 

7-8 

ended  confinement  systems. 

In  this  Letter,  we  present  measurements  of  the  ponderomotive  potential 

near  gyroresonance  and  compare  them  with  a  calculation  which  includes  non- 

adiabatic  effects.  It  is  found,  both  experimentally  and  theoretically,  that 

<J>p  is  vanishingly  small  at  w  =  fl  rather  than  singular.  In  fact,  the  maximum 

enhancement  in  $  does  not  occur  at  u  =  ft,  but  at  the  transition  from 
P 

the  nonadiabatic  to  the  adiabatic  regime,  namely,  for  w  ~  II  ♦  1.7v/£. 

The  experiment  is  performed  by  injecting  a  low  density  (single  particle 
limit)  monoenergetic  pulse  of  ions  through  an  electrostatic  rf  structure  as 
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shown  in  Fig.  1.  The  rf  electric  field  amplitude  is  increased  until  the  trans¬ 
mitted  current  disappears.  At  this  point  4^  is  equal  to  the  kinetic  energy 
2 

mv  /2  of  the  particles.  Experimentally  one  must  be  extremely  careful  that  the 
decrease  in  current  arises  because  the  particles  are  reflected  by  the  pondero- 
motive  force  rather  than  being  scraped  off  by  the  antenna  as  a  result  of  cyclo¬ 


tron  heating.  The  ion  source  is  a  thermionic  emitter  capable  of  producing  'Li 

27  + 

and  At  simultaneously.  A  pulse  of  ions  of  energy  q  Vg,  where  Vg  is  the  source 
bias  voltage,  is  produced  by  gating  a  grid  from  V  >  Vg  to  V  =  0  for  a  duration 
of  20y  s.  The  ions  are  guided  along  the  axis  of  a  grounded  conducting  cylinder 


by  a  uniform  axial  magnetic  field  Bq  =  1.8  Kg.  The  beam  radius  is  0.5  cm,  and 
the  cylinder  radius  is  7.3  cm.  The  rf  antenna  consists  of  two  electrically 


isolated  electrodes  made  by  splitting  a  cylinder  in  half  lengthwise.  It  is 
excited  through  a  center  tapped  transformer  so  that  E  is  perpendicular  to  B 
on  axis.  We  have  calculated  the  vacuum  rf  field  E(r,  e,  z)  analytically  to 


relate  E  to  the  rf  potential  +^Tf  applied  to  the  electrodes.  On  axis,  the 

*  2  2 

field  has  a  Gaussian  profile  E(r  =  0,  0,  z)  =  y  Eq  exp(-z  /2l  )  with  scale 
length  SL  which  can  be  varied  by  changing  the  length  of  the  antenna.  The  ions 
which  pass  through  the  antenna  are  collected  by  an  electrode  which  is  137  cm 


away  from  the  source. 

Figure  2  demonstrates  the  mass  selectivity  of  the  ponderomotive  force  as 

well  as  the  distinction  between  true  ponderomotive  stopping  and  fake  current 

drop  caused  by  resonance  heating.  The  time  dependence  of  the  collected  ion 

current  is  shown  for  various  conditions.  Since  all  emitted  ions  have  the 

same  energy  the  arrival  time  for  each  species  depends  on  the  mass  m.  Without 

7  + 

rf  the  signal  exhibits  two  peaks,  one  corresponding  to  in  *  7  (  Li  )  and  the 
27  + 

other  m  »  27  (  At  ).  When  the  rf  pulse  is  applied  at  the  gyrofrequency  of 
m  *  7,  i.e.,  u  ■  the  Li  peak  disappears  and  the  At  peaks  remain  unaffected 


(u)  >  ^27) •  The  drop  in  Li  current  is  not  a  ponderomotive  effect,  but  rather 
is  the  manifestation  of  ion  cyclotron  heating.  Our  measurements  indicate  that 
the  orbits  become  so  large  that  the  particles  are  scraped-off  by  the  antenna 
structure.  When  the  frequency  is  increased  to  u  s  1.5  as  seen  in  Fig.  2, 
the  early  peak  (Li)  disappears  once  again,  but  now  an  additional  peak  appears 
just  after  the  rf  is  turned  off.  These  ions  are  reflected  twice,  first  by  the 
ponderomotive  potential  at  the  antenna,  and  then  by  the  gate  grid  at  the  en¬ 
trance  end  of  the  device.  The  area  under  the  additional  peak  indicates  that 
90%  of  the  Li  ions  are  confined  by  the  rf.  The  A Z  peak  is  unaffected  because 
a)  >>  ^27  indicating  the  mass  selectivity  of  the  ponderomotive  force.  Similar 
mass  selection  occurs  near  fl  7. 

The  ponderomotive  potential  is  measured  for  a  fixed  value  of  w,  $1,  Z 
and  v  by  increasing  until  the  transmitted  ion  current  disappears.  For 
this  value  of  the  ponderomotive  potential  equals  the  initial  kinetic 
energy.  If  the  gyrofrequency  enhancement  of  the  ponderomotive  potential  is 
large  then  a  relatively  small  rf  electric  field  is  required  to  reflect  the 
particles  and  vice  versa.  Data  for  which  there  are  no  trapped  ions  are  dis¬ 
carded  because  the  previously  described  scrape-off  would  influence  the  measure¬ 
ment.  The  initial  kinetic  energy  of  the  ions  is  qVg  +  15%  as  verified  by 
time-of- flight  analysis  and  retarding  field  analyzer  measurements.  The  elec¬ 
tric  field  amplitude  on  axis  Eq  is  obtained  from  by  using  our  analytical 
solutions  for  the  vacuum  field.  This  is  justified  because  the  ion  density 
is  purposely  small  to  minimize  space  charge  effects.  The  gradient  scale 
length  Z  is  also  obtained  from  these  calculations  since  the  antenna  radius 
and  length  aTe  known. 


Figure  3  shows  the  behavior  of  the  ponderomotive  potential  near  the 

2  2  2 

gyrofrequency.  $  is  scaled  to  q  EQ/4mD  as  suggested  by  Eq.  (1)  to  fa- 
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cilitate  comparison  of  the  data  (dots)  with  the  theoretical  predictions  (solid 
curves).  The  curve  for  v/fcft  =  0  is  given  by  Eq.  (1)  and  the  rest  will  be  de¬ 
scribed  later.  Five  sets  of  data  for  different  values  of  l  and  v  are  shown 
to  demonstrate  the  important  parametric  dependences.  Away  from  gyroresonance 
4>p  agrees  quantitatively  with  the  prediction  of  Eq.  (1)  whereas  near  ft  the 

data  disagrees  even  qualitatively;  $  reaches  a  maximum  at  a  frequency 

P 

greater  than  ft  and  tends  toward  a  small  value  at  ft.  There  are  three  sets  of 
data  which  have  the  same  value  of  v/Jlft  but  different  values  of  i l  and  v.  These 
data  agree  with  one  another  within  the  experimental  error  indicating  that  the 
relevant  parameter  is  v/ift.  By  obtaining  data  over  a  wide  range  of  parameters, 
we  find  that  the  maximum  in  occurs  at  a  frequency  w  ;  ft  +  1.7  v/l.  As  a 
result,  the  maximum  enhancement  of  the  ponderomotive  potential  depends  on  var¬ 
ious  parameters  through  the  figure  of  merit  v/J Ift;  the  implication  is  that 
in  order  to  obtain  a  significant  enhancement  of  the  ponderomotive  force  near 
gyroresonance  v/Uft  must  be  small. 

To  explain  these  results  it  is  necessary  to  integrate  the  equation  of 
motion  properly  as  a  particle  traverses  an  electric  field  structure  of  finite 
extent.  Consider  the  exact  equation  of  motion  for  a  particle  with  position  r 
in  the  magnetic  field  B  =  Bqz 

r  *  E(r,t)  +  ftr  x  z  (2) 

_  m  _  _ 

A  A 

For  simplicity,  assume  that  E  =  (E  (y,z)y  +  E  (y,z)z)  cos  (u>t  +  i|>)  where  ♦ 

-  y  z 

is  a  phase  factor  and  V  *  E  =  0.  We  proceed  to  solve  Eq.  (2)  by  employing 

a  multiple  time  scale  perturbation  expansion  r(t)  =  rg(t)  +  r^(t)  where  rf 

gives  the  fast  response  of  the  particle  and  r  varies  on  a  slow  time  scale 

„s 

(compared  to  27r/w) .  By  expanding  E(r,  t)  about  r  assuming  |rj  «  i  and 
equating  terms  of  similar  time  scales  we  obtain 
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r,  =  ^  E  (r  ,  t)  +  Q  r.  k  z 
_f  m  _  \s  „f 


(3) 


where  E(rs,t)  is  evaluated  at  the  average  location  rg.  The  slow  response 
of  the  particle  is  determined  by  the  static  ponderomotive  force 


fP  =  <(!f  *  t}> 


(4) 


where  the  brackets  denote  an  average  over  the  fast-time  scale  2tt/uj.  The  par¬ 
allel  component  responsible  for  stopping  the  particle  has  two  contributions 
<qyf  Ez>  and  <q  ~  Ez>.  The  latter  gives  the  usual  unmagnetized  con¬ 
tribution*  which  is  not  of  interest  here  and  is  quite  small  in  our  experi¬ 
ment  because  E^  =  0  on  axis  (y  =  0)  where  the  particles  are  injected.  Using 

V  x  E  =  0  the  first  term  becomes  <q  y c  E  >.  To  evaluate  this  we  solve 

..  M  /f  dz  y 

for  y^  by  integrating  Eq.  (3) 
ft 


y  =  S_ 

yi  mSl 


dt'  E  (0,  z  (t'))sin  (fl(t  -  t'))cos(a)t'  +  ty) 
y  s 


(5) 


In  the  adiabatic  limit  £(ui  -  fl)/v>>l,  E  (0,  z  )  can  be  removed  from  the 

y  s 

integral,  y „  becomes  singular  and  the  standard  result,  i.e.,  Eq.  (1),  is 
recovered.  However,  near  gyroresonance  must  be  retained  inside  the  time 
history  integral.  In  general,  Eq.  (5)  is  diffifult  to  solve,  hence  we  make 

A 

the  simplification  that  z(t)  =  vt  z.  That  is,  we  integrate  over  the  unper¬ 
turbed  particle  orbit  retaining  E^  inside  the  integral.  We  shall  return 

to  this  point  later.  Substituting  y^.  into  Eq.  (4)  and  defining  the  pondero- 

[z 

motive  (pseudo)  potential  as  <j>p(o,z)  =  -  I  dz'  (o,z')  •  z  we  obtain 


^(0*2)  -  gmuJ^u> 


8muiv 


4mwv 


(6) 


where  Im  refers  to  the  imaginary  part  and 


In  the  adiabatic  limit  (o>  -  ft)  £/v  >>  1,  Eq.  (6)  reduces  to  Eq.  (1).  Otherwise 
the  strength  of  the  potential  has  an  additional  dependence  on  the  electric  field 
profile  and  the  velocity  of  the  particle. 


Using  a  Gaussian  profile  for  E^  in  Eq.  (6)  gives  an  analytical  expression 

Q 

for  <j>p(o,z)  in  terms  of  the  plasma  dispersion  function.  Proceeding  to  identify 

<!>p(0,0)  with  the  maximum  stopping  power  measured  in  the  experiment,  one  arrives 

at  the  solid  curves  exhibited  in  Fig.  3.  The  agreement  between  the  calculation 

and  the  experimental  values  is  excellent.  This  result  also  agrees  with  previous 
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rf  plugging  experiments  in  which  the  maximum  stopping  power  of  the  pondero- 
motive  force  occured  at  frequencies  above  the  gyrofrequency. 

It  is  illustrative  to  evaluate  <|>p(0,0)  for  the  Gaussian  profile  near 
gyroresonance,  i.e.,  (oj  -  ft)«./v  <<  1, 


2£2  2g2^2 


(8) 


which  is  finite.  The  first  term  is  small  and  arises  from  the  nonresonant 
polarization  of  the  antenna,  i.e.,  the  rf  field  is  linearly  polarized.  The 
second  term  arises  from  the  resonant  polarization  and  vanishes  exactly  at 
gyroresonance.  Physically  this  occurs  because  the  transit  time  t/v  is  very 
short  compared  to  the  time  required  for  the  particle  to  see  a  change  in  the 
relative  phase  of  the  perpendicular  rf  electric  field,  i.e.,  (id  -  ft)"*.  Con¬ 
sequently,  the  particle  response  depends  on  the  phase  of  the  electric  field 
when  it  enters  the  rf  structure.  In  this  case,  the  notion  of  a  time  indepen¬ 
dent  conservative  force  breaks  down  and  stochastic  particle  heating  occurs.*® 
This  effect  is  included  in  the  solution  of  Eq.  (3)  and  manifests  itself 
through  the  integral  £1  appearing  in  Eq.  (7) . 


It  is  somewhat  surprising  that  there  is  such  good  agreement  between 

theory  and  experiment  because  our  calculation,  as  well  as  the  adiabatic 
1-3 

calculation,  uses  the  Born  approximation.  This  assumes  that  the  high 
frequency  response  of  the  particle  can  be  calculated  by  integrating  the 
equation  of  motion  along  the  unperturbed  orbit;  a  condition  which  is  vio¬ 
lated  in  rf  plugging  experiments  including  ours  because  the  particles  are 
reflected  at  some  point.  To  further  validate  the  perturbation  results  we 
have  calculated  the  exact  particle  orbits  numerically  using  a  Gaussian  pro 
file  for  the  electric  field.  The  stopping  potential  is  also  found  to  be 
iiN^yse  agreement  with  Eq.  (6). 


In  summary,  a  second  order  perturbation  calculation  illustrates  the 
finite  nature  of  the  ponderomotive  effect  near  the  gyroresonance  and  is  in 
excellent  agreement  with  detailed  experimental  observations  with  test  ions 
Near  the  gyrofrequency  the  ponderomotive  potential  is  a  function  not  only 
of  position  but  also  depends  on  the  velocity  of  the  particles. 
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Figure  1  Schematic  of  the  experiment. 

Figure  2  Collected  ion  current  versus  delay  time  for  three  conditions 

of  the  rf  field.  Ion  energy  is  10  ev. 

Figure  3  Scaled  ponderomotive  potential  versus  frequency  for  different 

conditions.  Dots  are  measured  and  the  curves  are  the 
theoretical  predictions.  The  curve  for  v/£,ft  =  0  is  given  by 
Eq.  (1). 
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